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1. Introduction 

In the most general setting one can define a Gray code as a listing of words, formed 
with letters from an alphabet with k symbols, ordered in such a way that any two 
successive entries differ in just one symbol. These codes are named after Frank Gray, 
an engineer at Bell Laboratories who employed it in a patent in 1953 PO], although the 
use of the code has been traced back at least to 1878 when the French engineer Jean 
Maurice Emile Baudot (after whom the name baud for the signaling rate measure was 
chosen) developed a telegraph code on these grounds which readily prevailed in France. 

We are interested in the case where numerical symbols are involved. We can 
still distinguish different types of codes according to the number system used. The 
most commonly employed basis are of course 2 and 10 giving rise to the binary and 
decimal Gray codes respectively. Originally, the characteristic feature of these codes 
was important in mechanical devices such as counters, for no more than one piece 
at a time needs to act in every counting operation. Needless to say, this feature 
decreases the possibility of occurrence of errors albeit these may still be large. Although 
nowadays perhaps old-fashioned, a very illustrative and familiar situation comes from 
the mechanical odometer of a car. When it reads 9999 km. five inner wheels have to 
rotate to show 10000 km. Gray code is designed to avoid such multiple changes. 

In a wide variety of types Gray codes have been and still are applied in a large range 
of extremely distinct fields. We leave apart their use as puzzle-solvers |7j (e.g. the Towers 
of Hanoi or the Chinese rings) or in some even more bizarre application like campanology 

As a sample of more canonical fields where Gray codes are relevant we can cite as 
mere examples: database searching [T7], design of charged particle detectors ^|, space 
filling Hilbert curves genetic algorithms ^2] and their use in parameter optimization 
problems [211, neural networks [121, computation of unstable periodic orbits in chaotic 
attractors [Ij and, very recently, quantum computing [2S|- Even a little surfing in the 
web will convince the reader of their importance in commercial matters related to A-D 
converters or design of communication codes. 

In our case it was using the native Gray code to generate Walsh functions that we 
were led to the considerations presented in this paper. We take a Gray code as the 
play ground in which we build and represent some construction that is then analyzed 
by techniques familiar in other fields like dynamical systems or time series analysis but 
whose application is novel, so we think, in this context. The construction we refer to is a 
permutation, Pq, of natural numbers which in fact constitutes a known sequence, listed 
as A003188 in OEIS [23]. Our approach is basically descriptive and is not directly related 
to the interesting topic of the mutual relationship between Gray codes and permutations 
which has been extensively considered in the literature ^3] . In Section |21 the specific 
type of Gray code which we will focus our attention on is briefly reviewed, namely the 
so-called binary reflected Gray code (BRGC for short). Then we build what we call Gray 
curve which we analyze from a fractal geometry point of view in Section [HI Statistical 
considerations borrowed from modern time series analysis are applied in Section |3J For 
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easier reading we have postponed to Section El the discussion of the technical details of 
the structure of the permutation of 2^ elements associated to the native Gray code. In 
the last Section we add some final comments. 

2. The binary reflected Gray code 

The general definition of Gray code recalled in the Introduction corresponds to what 
is known in the literature as a (A^, /c)-Gray code and can in fact be fulfilled by very 
different orderings. Consequently one can find in the literature a considerable variety of 
members in the family of Gray codes and generalizations [221 ■ only interested 

in the length-iV binary subclass which includes all 2^ words of bits ordered in such 
a way that any two consecutive words differ in just one bit. We will even restrict our 
attention to the binary reflected Gray codes which are formally defined in a recursive 
way ^ : 

Let Tq denote the empty string and F^r be the BRGC sequence of A^-bit strings, 

then 

Tn+1 = oTn, ir^ (1) 

where OTn denotes the sequence Fat with prefixed to each string, and IF^ denotes 
the sequence F^v in reverse order with 1 prefixed to each string. 
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Table 1. Native Gray code construction for the first 2^ natural numbers and its 
associated permutation Pq. 

Perhaps the easiest way to seize the idea that underlies the reflected binary Gray 
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code is through an example. Table-Q shows the first 16 decimal integers 0, 1, . . . , 15 and 
their corresponding binary Gray code with = 4. The obvious first column lists the 
natural decimal ordering we are familiar with. The second column, which we shall not 
use explicitly in the following, is included just to explicitly appreciate the differences 
between standard base-2 and Gray code orderings. The latter is shown in the fourth 
column. A first look at it reveals the aforementioned defining feature: adjacent entries 
differ in one, and only one, bit. The third column is incorporated to illustrate more 
clearly the building procedure from equation (0). The fifth column, which is central to 
our discussion, displays the decimal version of the previous column. 

We have identified four blocks of rows in the Table-Q which will help in the 
explanation that follows, referred to the third column. The first block with two rows 
can be considered the seed for the Gray code, it represents Fi from equation (jT]). The 
next block includes the previous one together with the added next two rows. When 
observing only the right-most bit in these latter one sees simply the previous two rows 
in reverse order, i.e. Ff . When and 1 are added respectively to the left we complete 
F2. Observe that in order not to obscure the illustration of the building process the 
has not actually been written in the first two rows. The third block adds 4 new rows 
which accommodate the first 4 rows in reverse order and then the left-most 1 is prefixed. 
The fourth block illustrates the procedure again to obtain F4 which is given in full in the 
fourth column. The generalization to 2^ elements is straightforward. This is a recursive 
algorithmic way to construct the Gray code table. 

There are other possibilities for describing an algorithm leading to Gray codes. For 
example, once we have chosen the number of bits, A^, we can start with an arbitrary 
string of length A^. The next term in the sequence is obtained by changing only one 
bit. From here on the same procedure is repeated with the proviso that strings which 
have already appeared in the listing are discarded. Observe however that in this way we 
obtain certainly a binary Gray code, although we cannot assure it is the native BRGC. 

Apart from these recursive and constructive ways to obtain a BRGC of given length 
there is also an explicit method Let us write the kth entry of the BRGC of length 
A^ as the string 6^(/c)6^ (/c) . . .b^{k), then the jth bit (from left to right) is given by 
the parity of a binomial coefficient by the following expression which will be helpful in 
Section-|ni : 

bfik) =2--2---i q2._,._,_,„(,^,)/2] (mod 2) (2) 

where [■] is the integer part function and l<j<A^, 0<A;<2^ — 1. A comment may 
be in order here. The previous equation could give the wrong impression that the Gray 
code representative of the natural number k may sensibly change with the total number 
of bits A^. Either from equation or directly from one can easily see that this is 
not the case. Increasing A^ only adds O's to the left of the representative string. 

The native code is cyclic in the sense that 16, in our example, would go onto 0. 
More generally, the code is cyclic provided we deal with a number of elements that is a 
power of 2. It is worth mentioning here that the exact number of cyclic Gray codes for 
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general A^, which is related to the number of hamiltonian paths in an hypercube, is still 
an open problem in combinatoric Analysis. From the numeric stand point, the Gray 
code listings that we use here were generated with the short Fortran routine in pDl . 

Eventually, we are particularly interested in interpreting the entries in the Gray 
code column as natural numbers written in base-2. Their decimal conversion is given 
in the last column of the Table-^ Comparison of first and fifth columns admits the 
following interpretation which we word in two obviously equivalent ways: 

(i) The last column in Table-Q is a specific permutation of the 16 elements 
{0, 1, 2, ... , 15}. Or in general, the length N BRGC gives a permutation of the 
2^ elements {0, 1, 2, . . . , 2^ — 1} which we denote by Pq- 

(ii) The last column in Table-Q defines a function of a discrete variable in a definition 
domain with 16 (in general 2^) points. 

The reason we have stated separately these, otherwise interchangeable, sentences is 
because each of them suggests a different line of approach. In fact, while the permutation 
terminology is more amenable to combinatorial analysis, the functional sense calls to 
mind the idea of discrete dynamical system associated to the iteration of functions and 
related fractal structures. Both interpretations will be considered in this paper and 
yield interesting findings when considering 2^ elements with N ^ 1. By now, we 
stress the fact that all of our development pivots about the correspondence between the 
first (decimal sequence in standard order) and the fourth columns (decimal equivalent 
of the binary Gray code entry) in Table-^ Gray codification has been then just an 
intermediary in the process. 

We introduce now some notation. Let = = n, n = 0, 1, 2, 3, . . . , 2^ — 1} 
denote the sequence of the first 2^ natural numbers in base-10 and 0^^' = {(?„,, n = 
0, 1, . . . , 2^ — 1}, the decimal representation of the Gray code, i.e., for = 4, the last 
column in Table-Q Pq denotes the permutation that connects 2}'^^ and 0'^^ for a given 
A^, i.e., Pg is the element of the symmetric group which acts on S!)'^^ according to 

PGin) = Gn, n = 0,1,..., 2^-1. (3) 

In the next two Sections we present a series of heuristic considerations on Pg deferring 
to Section El its more technical details. 

3. A self-similar Gray set 

Let us now represent graphically the action of Pg- We consider then the set SIat 
of the plane consisting of the points with coordinates (n, PG'(n)) = (-D„,,G„), 
= 0, 1, . . . , 2^ — 1. In Figure-^ we plot ©'^^ versus S}'^^ up to A^ = 18, where we have 
connected successive points only as a visual help. A direct view of the disconnected set 
is given in the lower inset of the same figure. When no confusion may arise we shall 
continue to refer to this set as the Gray code curve. We could describe it, at first, as a 
kind of crinkly devil's staircase. In view of the algorithmic way explained above to build 
up iSit^l, it becomes clear the self-similar nature of the plot: Increasing a new bit to the 
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leftmost position of the Gray number means to sum a certain power of two, and besides, 
the preceding binary structure is inherited with a reversed ordering. In the log2 — log2 
representation of the upper inset the scale invariance is fully appreciated. In the present 
section we describe the set 'Qi^ for large from the point of view of fractal geometry. 

3.1. The Gray code curve as a fractal 

The self-similarity of the Gray code curve may be made explicit. Figure-|21 contains its 
first 32 points. Two subsets of 16 points each have been framed. The subset of points 
in frame B is geometrically obtained from subset A by a procedure of the type replicate, 
reflect and translate, as follows: 

(i) Make an adjacent copy of the set of points in panel A. 

(ii) Reverse the abscissae order in the replicated copy, i.e. substitute it by its image in 
a vertical mirror. 

(iii) Add 2^ (= cardinal of the subset) units to every ordinate. This is panel B. 

It is easy to see that inside plot frame A a similar procedure has been followed starting 
with the eight first points. The seed of the procedure was just 2lo = {(0, 0)}. 
To put this description in formulae let 21^ be 

2tfe = {(n,PGH),n = 0,l,...,2^-l}c M^ (4) 

the set of the first 2^ points of our curve, k = 1,2, . . . Then 

=2lfcUFfe(2tfe) (5) 

where ^^(21^) is the image set of 21^ under the function 




(6) 



X \ _ / -1 \ I X 

y ) ~\ I )\y 

defined on which corresponds to the previously described manipulations. By this 
iterative procedure the whole curve in Figure-d grows out. From a purely geometric 
point of view the previous function defines an afiine transformation. But its translation 
part depends on k, the generation index, a feature to which we shall come back. 

Alternatively we can inscribe the curve in the unit square with vertices 
(0, 0), (1, 0), (1, 1), (0, 1) by the following iterative procedure. Let us denote by *Bo the 
set {(0,0)}, incidentally 5Bo = 2lo, and define 

^k+i = M(^fc) U Lki^k) (7) 

where M(Q5fc) and L^i^k) are respectively the image sets of 05^ under the affine (still 
fc-dependent in the case of Lk) transformations 

(8) 
(9) 
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which scale down, reflect and translate the curve. From all these considerations the 
fractal flavor of the whole construction is evident. 

3.2. IFS generation of the Gray curve 

The procedures just described to construct the Gray set are clearly reminiscent of 
the so called Iterated Functions Systems, very popular in fractal building and image 
compression . As it is well known a collection of functions are applied, either following 
a deterministic plan or randomly, to a simple object. The two more striking features 
of this approach are, on one hand, how astonishingly complex images arise from very 
simple seeds. On the other hand perfect coincidence in the final product of deterministic 
and random processes, although more and more familiar every time thanks in part to 
the popularization of chaos studies, is always odd. Although not necessarily so, in most 
of the applications of IFS the functions involved are a few affine transformations with 
constant parameters. 

It is natural then to look for an IFS generator of our set SItv- The problem 
is, however, that one cannot find constant parameters for a bidimensional affine 
transformation capable of generating by iteration our fractal. The reason being the 
k dependence mentioned earlier. So we are led to look for more general functions, i.e. 
non-linear ones, of the coordinates, such that when iterated produce the desired result. 
But the k dependence continues to be a nuisance. The only way out is to increase the 
dimension of the space we are working on from two to three. This is exactly the same 
idea one takes advantage of in, e.g.. Classical Mechanics to convert a time dependent 
hamiltonian system in an autonomous equivalent. 

We then propose to extend to and take as our IFS the following 
functions: 



The price one has to pay for getting rid of the k dependence is the higher dimension 
we work on and the non-linear character of the function g involved. Needless to say 
that for representation purposes we only plot the projection on the xy plane. Anyway, 
equation ((7j) changes now to 



Observe that the action of / and g on the z coordinate prevents the existence of fixed 
points when each of these functions is iterated on its own. In the xy plane when / is 
repeatedly iterated the results tend rapidly to the origin. It is only a little bit more 
complicated to see that iteration oi g leads in the limit to the point (1, 1). The combined 
and competing iteration of the two functions is what produces the final figure. The set 
{f,g} constitutes a bona fide IFS. The system is deterministic if / and g are repeatedly 
applied one after the other. When this is done starting with (0,0), then the whole set 




(10) 



(11) 




(12) 
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21 AT is generated. It is a random IFS if / and g are iterated with prescribed probability 
i on an arbitrary initial point in the unit square. 21 at is obtained as the final attractor 
of the process. 

4. A time series from BRGC 

Our proposal now is to interpret G„ as successive values in a time series. Conventional 
statistics (average values, standard deviations and other indicators calculated from 
higher moments) are of no particular interest for us here because they depend solely of 
the numerical values G„ and these are simply the first 2^ natural numbers. The order 
in which the different values appear does not play any role in the computation of those 
statistics. The information from the time series obtained in this way could be termed 
as static. Instead, we are interested in the dynamics represented by the series. The 
ordering, then, matters and we need other indicators. With this aim we analyze in the 
present Section the series of increments {A„} with A„ = Gn — Gn-i, n = 1, 2, . . . 

4.I. Statistics on An 

The inset in Figure-Ol shows a sample of the series of increments {A,„}. By its very 
definition the average of the first m values is Gm/^ and because it is always positive 
the Gray curve is globally increasing. The statistics on the set of values {|A„|} reveals 
that their frequency obeys a power law, with the exception of the very small values 
(see main graph in Figure-El). As a matter of fact, the frequency histogram goes like 
~ 1/1 A|, which means that the probability of encountering a jump |A| < 5 grows as 
log 5. One could express this fact by saying that the values Gn are far from equidistant. 
By contrast, it seems appropriate to mention here that in its binary form, BRGC is 
characterized by any two consecutive elements being equally spaced by a unit Hamming 
distance. 

4-2. Residence times from A„ 

An often encountered statistic when analyzing time series is the so-called residence 
time, or sojourn time, r. To deal with it, one has to define a dichotomic property of 
the system such as crossing a predefined threshold and keep track of the time elapsed 
between successive crossings. In our case we focus our attention on changes in sign of 
A„. Accordingly we define the residence times as the length of the sequences in which 
A„ keeps its sign. If it stays just one unit in the, say, positive region, we define its 
sojourn time to be just one unit. 

When we collect the statistic for residence times in the set of values of A„ the 
fact which first calls our attention is that r can only take three different values: 1,2 
or 3. Suppose A„ has changed sign with respect to A„_i, then there are only three 
possibilities: either A^+i has sign opposite to that of A„ (then r = 1), or A.„+i and A„ 
have equal sign but different from that of A„_|_2 (which corresponds to r = 2), or A„, 
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A„+i and A„_|_2 have the same sign which necessarily changes in A„+3 (giving r = 3). 
Never more than two consecutive steps can be given without encountering a change of 
sign. As a matter of fact the statistics reveals more than that: the probabilities (relative 
frequencies of appearance) for the values of r are exactly 

p{l) = 0.25,j9(2) = 0.5, p(3) = 0.25, p(r > 3) = 0. (13) 

This behavior of the residence times contrasts with other known situations. For example, 
exponentially decreasing probability distributions are common in chaotic systems. 

It is worth mentioning that the statistics of residence times we are treating provides 
information only about changes of sign in A„. It says nothing about the direction 
(upwards or downwards) of the change. The fact that, as mentioned earlier, the average 
of increments {A„} is positive suggests the increasing overall character of our curve. By 
simply looking at it at any scale one corroborates that it is globally increasing, in spite 
of its ups and downs we perceive in more detailed graphs. In the jargon of time series 
analysis one says that this reflects extremely persistent dynamics which we analyze next 
from the Hurst exponent point of view. 

4-3. Hurst's R/ S analysis of An 

The so-called re-scaled range statistical analysis, R/ S analysis, is a method designed to 
reveal long-run correlations or anti-correlations in a complex process. The idea was first 
introduced in the analysis of some hydrologic time series, but has been also applied in 
other studies in climatology, geophysics, econometrics,. . . 

Following Feder we define the average increment over a period of k units as 

k 

{^)k = lj2^r = Gk/k (14) 
r=l 

The local accumulated departure of A„ from the mean (A)^ is then 

n 

X{n, k) = Y}^r - (A)fc] = G„ - ^Gfc (15) 

r=l 

The range of X(n, k) over this period reads 

i?(A;) = maxX(n, A;) — minX(n, fc) (16) 

where the extremal functions are evaluated in the interval 1 < n < k. The range 
increases with k. When properly normalized to the standard deviation 

1/2 

^ = I r > :\^r - {A)uf 1 (17) 



then Hurst found that the re-scaled range R/ S is well described by the power law 

R/S = (A;/2)^ (18) 

where H is called Hurst exponent. Purely random series give H = 1/2. Curves with 
H > 1/2 {H < 1/2) witness a long run correlation of persistent (anti-persistent) 
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character. This feature means that positive changes are more hkely followed by positive 
(negative) ones, whereas for random phenomena both results for the next move are 
always equiprobable. Hurst found that most natural phenomena give rise to 
~ 0.7, which is sometimes referred to as Hurst phenomenon. 
In the case of the Gray curve a Hurst-like analysis provides a value H = 0.97(3), 
revealing the extreme persistence of the series of increments, as pointed out above. This 
analysis appears in Figure-IU where R/S (straight line with open circles) is to be read 
on the right scale. 

4.4- Fractal dimension of the Gray curve 

We follow a technique by Higuchi jTT] to measure the fractal dimension of a set of points 
{n, Gn}, n = 0,1,2, ... , 2^ , forming the graph of a function like that in Figure-^ The 
idea consists in measuring the average length {L{k)) of samples with k points each. 
Then, if {L{k)) oc k~'^, the curve is fractal with dimension d. For instance, a time series 
obtained with pure noise should give d = 1.5, whereas for a regular curve d = 1. The 
method recalls the so-called Richardson plot in fractal geometry. 

In Figure-Hl we have plotted in doubly logarithmic scale {L{k)) against k, with 
= 18 (straight line with solid dots). The straight line is fitted to the points by the 
least square method and gives a slope d = 1.09(1). Here the error estimate comes from 
the stability of the result with respect to the number of data rather than from the fitting 
procedure. 

These values for H and d are consistent with the connection formula relating H 
and d which reads: d = 2 — H. 

4.5. Return plot of the Gray curve 

A common tool in nonlinear analysis of time series consists in giving a 2D- representation 
of the ID time series called return plot. A return plot, or first-return map, of the 
sequence {Gn} is a 2D-plot of the points {((?„_!, (?„)}. Thus, for instance, a true 
random sequence gives rise to a uniformly scattered plot. Chaotic sequences where the 
generation of the nth term involves only the (n — l)th term produce ID-patterns. In 
the present case, the return plot for the Gray curve is in Figure-|H1 We have drawn not 
only the points but also a path connecting them in order to realize the variety of jumps 
involved, in accordance with the set of values {A„}. The inset in Figure-El is a zoom of 
the indicated zone showing the detail of the curly structure of the curve. It is clear how 
strongly the natural numbers sequence gets intermingled in the permutation Pq'. The 
recurrence plot for the set {-Dn} is simply a straight line like y = x + 1. 

4-6. Flying distances in the random IFS Gray fractal 

At the end of subsection 13.21 it was pointed out that the Gray fractal can be obtained as 
the attractor of a random process of iteration of functions / and g. Here we investigate 
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the statistics obeyed by the distances traveled in every move. This point of view quits 
the static perspective on the fractal and enters a new kind of dynamics. 

Starting from an arbitrary point in the unit square, runs of 2^^ points were carried 
out. We have retained just the last 2^^ points, and the previous ones were considered as 
a transient. The idea is to lose any memory about the initial position by getting close 
enough to the true attractor. A histogram of the Euclidean distances traveled between 
two consecutive points in a random walk on the Gray fractal is shown in Figure-IHl It is 
noteworthy the presence of equally spaced peaks on the log-scale. The same is true for 
the bottom of valleys. This phenomenon means that the geometry of the fractal selects 
(or avoids) some particular flying distances along its random generation. The specific 
mechanism underlying this selection rule is unclear, as according to the statistic of A„ 
in subsection- l4.1l all length distances are present in the fractal itself. 

5. The native Gray permutation 

In this section we collect some properties of the permutation Pq which have been 
instrumental in most of the considerations presented so far. In correspondence with 
the name native Gray code sometimes used for the BRGC we call Pg the native 
Gray permutation. A little experimentation and some reflection lead to the following 
observations: 

(i) The recurrent definition of the permutation Pq as a function with domain S}'^^ is 

Pg(2" + m) = 2" + Pg(2" - m - 1) (19) 

with -Pg(O) = 0, and where n = 0, 1, . . . , — 1, m = 0, 1, . . . , 2" — 1. Particularly 
interesting values are Pg(2" - 1) = 2"-^ and Pg(2") = 3 ■ 2"-\n > 1. 

(ii) An explicit expression for Pcin) can be deduced from equation Q. Given n G N, 
n > 1, the minimum number of bits necessary to write it in base-2 is v{n) given by 

v{n) = [\0g2in)] + 1. 

Then we have 

v(n) 

Pc(n) = ^6;("^(n)-2^(")-^' (20) 

where h"j^"\n) has been defined in Section-|2 

(iii) The equation Pg{x) = x has, for any value of A^, only two solutions: x = and 
x = 1. It is to say, and 1 are the only elements of St^' which Pg leaves invariant. 
All other elements of under the action of Pg are grouped in cycles, but only 
cycles of length 2' will appear. In symbols that means that we can factor out Pg 
in cycles in the form 

k{N) 

Pg=11 cT (21) 
1=0 
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where q represents a cycle of 2' elements of which ni different specimens appear in 
Pg- k{N) indicates that the length of the longest cycle or cycles is 2'^^^\ 

(iv) The distribution of the 2^ elements of S'^^ in cycles follows a strict procedure: 
they fill cycles q in increasing order of /. No 2'— cycle appears until all allowable 
2'"^— cycles are filled up. Moreover, rti has a maximum value independent of A^. In 
other words, once a certain number of 2'— cycles has been reached no matter how 
big N becomes no other 2'— cycle will be formed. Furthermore, when N increases 
the factorization of Pq keeps the structure corresponding to the lower values of A^. 
The new cycles formed are as long as the longest ones for previous values of or 
longer. This property confers to Pq a certain type of shell structure. 

(v) 2'— cycles will appear for the first time for = 2'^^ + 1. In the language of the 
theory of dynamical systems one could say that at these values of A^ a period 
doubling bifurcation appears. This way of speaking has to be considered simply as 
an analogy and cannot be taken too literally. It is illustrated with an example in 
Figure-|Z| which is self-explanatory. 

With all these considerations in mind a careful book-keeping shows that the length 
of the longest cycle in Pq is 2^^^^^ with 

'\n{N -I) 



k{N) 



In 2 



+ 1 = [log2(Ar -!)] + !, N>2 (22) 



where [■] stands for integer part. Obviously for A^ = 1 we have k{l) = and, for any A^, 
no = 2. 

According to observation (iv) above, the maximum value of rii, n™"^ is, for / < k{N), 
independent of A^ and it is not difficult to see that 

^max _ 2(2'-'-0(2(2'-^) _ 1), Q < I < k{N) (23) 

while the number of longest cycles is 

p =N-2^-^ (24) 

where in these formulae k{N) has been for simplicity written as k. Observe that in this 
last case the dependence on A^ is explicit. As a test one can check that 

k{N) 

2^ni = 2^ (25) 

1=0 

To fully appreciate the complexity of the cycle structure of the Gray permutation 
Pg we show in Table-El a few cases of values characterizing that structure. It gives the 
complete cycle structure of Pq up to A^ = 10 and includes also the total number of 
cycles, Nc, and the average cycle length, lav 

The reminiscence of the atomic shell-model mentioned earlier is of course purely 
symbolic but it is reinforced when observing Figure-|H1 In it we represent k{N) as a 
function of A^ and one can see that the jumps occur for values of A^ which are exact 
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6.4 
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7.11 
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16 


52 


9.85 


10 


1024 


2 


1 
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30 


48 


84 


12.19 



Table 2. Cycle structure of permutation Pq 

powers of 2. The inset shows that precisely at these values the number of longest cycles 
falls down. 

Of course the general formulae we have deduced do not give explicitly the elements 
in each cycle like for example, for = 4, 

Pg = (0)(1)(2, 3)(4, 6, 5, 7)(8, 12, 10, 15)(9, 13, 11, 14), 

but this information is much less relevant for our purposes. 

It is interesting to notice the somewhat distinguished role of Pq in the symmetric 
group S2N. It is clear that from the information given in this section one can easily 
draw the Young diagram corresponding to the equivalence class containing Pq in the 
group S2N. It is also obvious that binary Gray codes other than BRGC originate also 
permutations in S2N, but these may have a cycle structure very different to that of 
Pg and consequently behave in other ways. On the other side any element Q G 5'2jv 
generates trivially a binary code, but in general it is not of the Gray type. 

6. Final comments 

Among combinatorics problems with physical applications the study of permutations 
has always been of particular interest. As recent exponents of it, in jB] and ^Hj random 
permutations are used in connection with stochastic processes. On the other hand Gray 
codes have also met applications in different branches of Physics, as widely shown in 
the Introduction. 

From the point of view of combinatorics or graph theory relationships between 
this type of codes and permutations are manifold and are being explored in different 
directions in mathematical literature [SlElllSl- But from physics standpoint no such 
kinship seems to have been studied in detail. In this paper we have analyzed the 
permutation associated with the binary reflected Gray code from the perspective, and 
with the techniques, of a dynamical system and its corresponding time series. In so 
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doing we have been led in a natural way to enhance some fractal aspects of the system. 
One may hope that this line of approach could be helpful in understanding a more 
dynamical role of permutations. 

The results of the analysis carried out in this paper involve only integer arithmetic, 
some book-keeping and graphical considerations and so constitute, in our opinion, 
appropriate material for pedagogical purposes. 

As a final summary we would like to mention as a particularly appealing result of our 
approach the generation of the permutation Pq by an iterated function system apphed 
either in a deterministic or a random way. Also worth mentioning is the discovery of 
strong selection rules acting on the system which manifest themselves in the distribution 
of residence times with only three allowed values, and of flying distances. Which, if any, 
of these results apply to other permutations associated with other codes may perhaps 
be worth considering in the future. As would be a deeper understanding of the selection 
rules mentioned before. 
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Figure 1. Self-similar Gray code curve 




Figure 2. Self-similar construction of the Gray code curve 
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Figure 4. Determination of the fractal dimension (left scale, solid dots) and Hurst 
exponent (right scale, circles) 
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Figure 5. 2D trajectory for the Gray curve and zoom of the marked area 
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Figure 6. Histogram of Flying Distances in the random IFS Gray fractal with 
and bin size 0.02 





Figure 7. Orbits for Gray permutation of elements to 31 




Figure 8. Statistic of cycles as a function of iV 



